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Warping and precession in galactic and extragalactic accretion disks 

Anderson Caproni^'^, Mario Livio^, Zulema Abraham^ and Herman J. Mosquera Cuesta ^''^ 

ABSTRACT 

The Bardeen-Petterson general relativistic effect has been suggested as the mech- 
anism responsible for precession in some accretion disk systems. Here we examine 
separately four mechanisms (tidally-induced, irradiation-induced, magnetically-induced 
and Bardeen-Petterson-induced) that can lead to warping and precession. We use a 
sample of eight X-ray binaries and four Active Galactic Nuclei (AGNs) that present 
signatures of warping and/or precession in their accretion disks to explore the viabil- 
ity of the different mechanisms. For the X-ray binaries SMCX-1 and 4U 1907-1-09 all 
four mechanisms provide precession periods compatible with those observed, while for 
CygX-1 and the active galaxies Arp 102B and NGC 1068, only two mechanisms are in 
agreement with the observations. The irradiation-driven instability seems incapable of 
producing the inferred precession of the active galaxies in our sample, and the tidally- 
induced precession can probably be ruled out in the case of Arp 102B. Perhaps the best 
case for a Bardeen-Petterson precession can be achieved for NGC 1068. Our results 
show that given the many observational uncertainties that still exist, it is extremely 
difficult to confirm unambiguously that the Bardeen-Petterson effect has been observed 
in any of the other sources of our sample. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — 
galaxies: nuclei — magnetic fields — stars: neutron — X-rays: binaries 

1. Introduction 

The importance of accretion processes in the line formation and continuum emission in AGNs 
and galactic sources, such as microquasars and X-ray binaries, has long been recognized in the 
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literature, even though the presence of an accretion disk in some of those systems is inferred only 
indirectly. 

The continuous improvement in the capability of telescopes and detectors to obtain high- 
resolution images, and spectra with higher sensitivity, allowed for the direct probing of the physical 
characteristics of some accretion disk systems (e.g., Ray et al. 1996; Jones et al. 2000; Jaffe et al. 

2004; Fathi ct al. 2006). 

In contrast to the standard picture of a flat disk surrounding the accreting object, some galactic 
and extragalactic sources present signatures of warping and precession in their disks and/or jets 
(e.g., Margon 1984; Sillanpaa et al. 1988; Abraham 2000; Gallimore et al. 2004; Caproni, Abraham 
& Mosquera Cuesta 2006). The apparent lack of correlation between the orientation of the radio 
jets and the plane of the host galaxy's disk might also be attributed to warped disks (e.g., Schmitt 
et al. 2002). 

The Bardeen-Petterson effect (Bardeen & Petterson 1975), a gravitational perturbation pre- 
dicted by general relativity, is one of the physical mechanisms that has been proposed to explain 
warping and precession in accretion disks (e.g., Scheuer & Fcilcr 1996; Nelson & Papaloizou 2000; 
Fragile &; Anninos 2005; King ct al. 2005; Caproni, Abraham & Mosquera Cuesta 2006). In order 
to explore whether this mechanism is universal, in the sense that it can be acting in both galactic 
and extragalactic accretion disks, we have selected eight X-ray binaries and four AGNs that present 
signatures of warping/precession in their disks. We have analyzed individually the predictions of 
three additional precession mechanisms for the same sample of sources, in order to compare these 
predictions with those from the Bardccn-Pctterson effect. In Section 2, wc describe the accretion 
disk model used throughout the paper, as well as constraints on the basic parameters of the accre- 
tion disk. In Section 3, we present the four warping/precession mechanisms studied in this work. 
A brief introduction of the sample of sources, as well as their basic parameters (inferred observa- 
tionally) is given in Section 4. A comparison between the theoretical predictions and observations 
is carried out in Section 5, and conclusions follow. 

2. The accretion disk model 

We consider an accretion disk through which mass is accreted onto a compact object (a black 
hole or a neutron star) at a mass accretion rate M. The angular momentum of the accretion disk 
per unit area can be written as: 

Lair) = S(r)^2(r)r^ (1) 
where r is the radial distance from the compact object, S is the surface density of the accretion 
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disk (integrated over a scale height H), and $7 is the (relativistic) Keplerian angular velocity of the 
disk, given by (e.g., Abramowicz, Jaroszyhski & Sikora 1978): 



Kg 



(2) 



where a* is the ratio between the actual angular momentum of the compact object and its maximum 
possible value and i?g = GMj(? is the gravitational radius of the compact object, (where G and c 
are the gravitational constant and the speed of light respectively). Note that for r 3> -Rg or a* ~ 0, 
we recover the Newtonian angular velocity. 

If the accretion disk is not self-gravitating, the scale height of the disk can be determined from 
hydrostatic equilibrium to be: 



air) 



(3) 



where Cg is the sound speed, given by (e.g., Abramowicz et al. 1988): 



5 din Q{r) ui{r)Q.{r) 
3 dlnr a 



(4) 



where a is the dimensionless viscosity parameter introduced by Shakura &; Sunyaev (1973). The 
azimuthal kinematic viscosity of the disk, ui, is calculated from (e.g., Krolik 1998): 



1^1 (r) = - 



M 



dlnO(r) 
27rS(r) [ dlnr 



R„ 



n{R^ 



(5) 



where R^s is the radius of the innermost marginally stable orbit, assumed to be the inner radius 
of the disk^. 

In this work, we assume a power-law disk (e.g., Ostriker, Shu, Sz Adams 1992; Maloney, 
Begelman Sz Nowak 1998; Caproni, Abraham Sz Mosquera Cuesta 2004; Raymond, Quinn Sz Lunine 
2005; Caproni, Abraham & Mosquera Cuesta 2006): 



^For neutron stars with strong magnetic fields, the accretion disk might be truncated at a larger radius than -Rms, 
at the uiaguetospheric radius. Even though this might be the case for some sources in our sample, it should not 
change our results in any substantial way. 
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S(0 = ^o{-^ 

its 



(6) 



where Sq and s are constants. 

For accretion disk models available in the literature, —2 < s < 2, while Eq should be determined 
from some reasonable assumptions concerning each system. 

An extreme upper limit for So can be found by imposing that the accretion disk mass is 
lower than the mass of the compact object, such that 



< 



27rG 



R: 



s+l 



(rRout 
Rms 



(7) 



In order to estimate the limit on Eq using equation (7) , it is necessary to know the outer radius 
of the disk Rout, which is generally not well constrained by observations, especially for AGNs. For 
our AGN sample, wc have adopted i?out ~ 10^i?g based on the model of CoUin-Souffrin & Dumont 
(1990), except for NGC 1068 in which intcrfcromctric maser observations indicate that its outer 
radius is at about 1.1 pc (Gallimore et al. 2004). In the case of binary systems in which accretion 
occurs via Roche lobe overflow, we take i?out ~ 0.88i?L (Papaloizou Sz Pringle 1977), where the 
Roche-lobe radius Ri, is given by (Eggleton 1983): 



Rl = 



0.49g2/3 



0.6g2/3 + ln(l + gV3) 



Rt) 



(8) 



where Rps is the binary separation. 

Another estimate for Sq can be derived from the assumption that the accretion disk is self- 
gravitationally stable. This imphes that the Toomre parameter Q = Csfi/(7rGS) (Toomre 1964) 
must be greater than unity, which leads to: 



M 



1/3 



R. 



out 



-(s+3/2) 



71/3 

Q ' 



(9) 



with the dimensionless function Fq given by: 



Fq{a*, Rout) 



1 - {RnjRoui f mn..)/mout) 

1 + a* {Rout/Rg)~^^^ 



(10) 
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We can see that Fq ~ 1 for ii!out ^ Rms-, which is true for all the sources in our sample. The 
upper limit on Eq corresponds to the minimum between S^'' and S^. 

On the other hand, a lower limit for Sq can be established from the requirement that the radial 
inflow is (highly) subsonic. This gives: 



^0 



> 



1 3aM fR 



T2 lOTTCi?, V R 



1 + a* 



Rr, 



Re 



-3/2' 



(11) 



where the Mach number T = Wr/cg, and is radial velocity of the disk material. Although the 

requirement of a subsonic accretion inflow only implies that T < 1, current accretion disk models 
usually give T < 0.01 (e.g., Shakura & Sunyaev 1973; Abramowicz et al. 1988; Narayan &; Yi 1995). 
Consequently, we have assumed T = 0.01 in this work. 



We therefore take for the allowed range of Sq in this work: Sq" < Sq < min 



^0 ' 



We 



should note that for the irradiation-driven instability, another lower limit on Sq can be derived 
from considerations of the disk opacity (see § 3.2). However, in most of the cases analyzed in this 
work, Sq' has provided a more restrictive lower limit for Eq. 



3. Warp/Precession mechanisms 

There are four main mechanisms that have been suggested for driving warping and precession 
in accretion disks. We consider all of these in turn. 



3.1. Tidal forces of a companion object in a binary system 

The precession of an accretion disk can be tidally induced by a companion in a binary system 
(e.g., Katz 1973; Sillanpaa et al. 1988; Katz 1997; Romero et al. 2000; Caproni & Abraham 2004a). 

We consider a binary system with masses Mp and Mg (for the primary and secondary, respec- 
tively), separated by a distance R^. Prom Kepler's third law. 



Ri 



G(Mp + Mg) ^2 



47r2 



orb' 



(12) 



where Porb is the orbital period. 

If the orbit is non-coplanar with the accretion disk, torques are induced in the inner parts of the 
disk, producing precession. Taking the outer radius of the precessing disk to be -Rprec (^prec < -Rout; 
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Romero et al. 2000), the precession period, Pprec, is given by (Papaloizou & Terquem 1995; Larwood 
1997): 



where n is the polytropic index of the gas (e.g., n = 3/2 for a non-relativistic gas and n = 3 for the 
relativistic case), q = Mg/Mp and 6 is the inchnation angle of the orbit with respect to the disk 
plane. The negative sign in equation (13) indicates that the induced precession is retrograde, in 
the sense of being contrary to the rotation of the accretion disk. 

Combining equations (12) and (13) (with the condition -Rprec < -^^out) we obtain: 



where f{q) is the function multiplying Rpg on the right-hand side of equation (8). Note that the 
ratio between the precession and orbital periods depends (apart from 9) primarily on the mass ratio 
of the binary system. 



Radiation produced by the accreting compact object can modify the dynamics of an opti- 
cally thick accretion disk. Petterson (1977) showed that if the disk is warped and optically thick, 
the radiation pressure will produce nonaxisymmetric torques that will change the initial warped 
configuration of the disk. Pringle (1996) and subsequently Maloney, Begelman & Pringle (1996) 
showed that even for an initially flat accretion disk, the radiation torques can warp and twist it. 
Pringle (1997) extended those calculations to include self-shadowing effects due to the warps. This 
radiation-driven instability has been invoked to explain anomalies in the morphology and variabil- 
ity in a huge variety of astrophysical sources (e.g., Cliffe et al. 1995; Maloney, Begelman & Pringle 
1996; Southweh, Livio & Pringle 1997; Livio & Pringle 1997; Armitage & Pringle 1997; Wijers & 
Pringle 1999; Ogilvie & Dubus 2001). 

Maloney, Begelman &; Nowak (1998) generalized the isothermal disk model assumed in pre- 
vious papers, and considered power-law surface density distributions ranging from the isothermal 
case {s = —3/2) to radiation-pressure dominated disks {s = 3/2). Solving numerically the twist 
differential equation governing such disks (without including self-shadowing), they found that the 
critical radius at which the disk becomes unstable is: 




>3 



(13) 




(14) 



3.2. Radiation-driven instability 
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-^cr — X^j. 



2 cn^'^ 



) Rg: (15) 



where ij = V2lv\ is the ratio between the vertical and the azimuthal viscosities in the disk, e is the 
accretion efficiency and x^t is equal to 27r for s = —3/2, increasing monotonically to ~ 4.89l7r for 
s = 3/2. 

In the steady-state regime, the precession period, -Pprec> due to the radiation-instability is 
(Maloney, Begelman &; Nowak 1998): 



^prfc = (-) , (16) 



jrad 

(To 



where 8.74 x 10 ^ < a < 1 (lower and upper limits referring to s = 3/2 and -3/2 respectively), and 
the dimensionless parameter ctq is given by: 



= 12c2Eti?|' ^^^^ 

where Lboi is the bolometric luminosity produced by the accretion process onto the compact object. 

As we have noted in Section 2, the disk has to be optically thick in order to be unstable to 
warping by irradiation. This implies the existence of an independent lower limit for the surface 
density, represented by Sq, which is not necessarily more restrictive than 

Note also that depends on which mechanism is the main contributor to the disk opacity; if 
Thomson scattering on free electrons of cross-section ctt dominates over free-free absorption (with 
opacity Kg), the optical depth can be calculated from tt = \/aTKsJrnE.Ti, otherwise rff = KffS 
(Shakura k Sunyaev 1973). 

The Rosseland mean opacity for free-free absorption is ks = koP^'"^/^, where p and T are the 
mass density and temperature of the accretion disk respectively, while is a numerical constant 
equal to 3 x 10^^ cm^ g~^ (e.g., Maloney, Begelman & Nowak 1998). Expressing p and T in 
terms of S and Cg, and assuming a non-relativistic ideal equation-of-state for the gas in the disk, 
Kff is given as: 



'5 A;b y^ nrmr) 
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where kB is Boltzmann's constant, fi is the mean molecular weight and mn is the atomic Hydrogen 
mass. 

If Kg < (TT/mu at the critical radius, the lower limit on Sq (in order to have an optically thick 
disk) is: 



where 



(20) 



and 



G-r (a* , Rc 



1^ 1^-'"'^'^" 
1 + a* I — 

^g 



FQ'/'(a*,i?cr). (21) 



On the other hand, if free-free absorption dominates: 

P \ -(^+3/4) 



(22) 



Thus, the lower limit for Eq in the case of irradiation torques is obtained from the more 
restrictive value between 'E'q (equation 11) and Sq. 



3.3. Magnetically-driven instability 

The influence of magnetic fields on accretion disks has been studied by several authors, espe- 
cially in the cases in which the central object is magnetized (e.g., Aly 1980; Lipunov & Shakura 
1980; Lai 1999; Terquem &; Papaloizou 2000). The production of quasi-periodic oscillations and 
luminosity variability in such objects has also been explored in the framework of external magnetic 
fields and accretion disk interactions (e.g., Agapitou, Papaloizou &; Terquem 1997; Terquem & 
Papaloizou 2000; Shirakawa & Lai 2002a,b). 
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Lipunov & Shakura (1980) and Lai (1999) showed that when the spin axis of the magnetized 
object, as well as its dipole moment, are not aligned with the angular momentum of the disk, 
magnetic torques will be generated, causing warping and disk precession. Extending those results 
to non-niagnctized central sources, Lai (2003) showed that a large-scale magnetic field (associated 
with magnetically driven outflows) threading the disk also induces a warping instability and a 
retrograde precession. Nonlinear simulations of warped, viscous accretion disks driven by magnetic 
torques were performed by Pfeiffer & Lai (2004). These numerical simulations showed that the disk 
can keep a steady-state warped shape with a rigid-body precession. 

In order to examine the applicability of the magnetically driven mechanism to our sample of 
objects, we have considered a very idealized configuration, similar to that adopted by Lai (2003). 
Let us assume that the accretion disk is threaded by a poloidal magnetic field Bp with a radial 
pitch angle if = arctan \ BYi/Bz\, where and Bz are respectively the radial and parallel magnetic 
field components in relation to the direction of the compact object's rotation axis (Z-direction). 

The poloidal magnetic field lines will be twisted by the disk rotation, generating a toroidal 
component B^ that has different signs above and below the disk plane due to the discontinuity of 
5r, at the disk mid-plane (Lai 2003). Following Lai (2003), we introduce the azimuthal pitch (, 
such that 5$ = tC-^z, with the negative and positive signs referring to B^ above and below the 
disk mid-plane, respectively. 

In addition, the accretion disk plane does not have to be perpendicular to Bz, being instead 
tilted by an angle (3 (see figure 2 in Lai 2003). In this case, the projection of Bz onto the perpen- 
dicular direction (to the disk) will be Bz cos (3, while the toroidal magnetic field will be ^(Bz cos (3. 

The discontinuities in the radial and toroidal components of the magnetic field inside the disk 
produce a net disk surface current. The interaction between the current's toroidal component and 
the radial component of the magnetic field leads to the magnetically driven precession, while the 
interaction between the toroidal component of the magnetic field and the radial component of the 
net disk surface current warps the disk. The magnetic torques per unit area responsible for the 
disk precession and warping (averaged over the azimuthal angle 0) were calculated by Lai (2003): 

{T^r^J)<t> = -'-^rBl (23) 



{T^^lU = -^rBlcosf3. (24) 

We can see that this mechanism produces retrograde precession, as in the case of the tidal 
torques in binary systems, and it pulls the normal to the disk plane away from the Sz-direction, 
increasing the angle P with time. 
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The precession period induced by the magnetic torques can be estimated by: 



where Jd is the integrated angular momentum of the accretion disk: 



"-Rprcc 

Jd(i?ms,i?prec) = 27r / La{r)rdr. (26) 



Rms 



The magnetic warping timescale can be calculated using: 



pmag _ JdjRms, -Rprcc) ,2'7\ 

warp jn^^^^^T^^.^^rdr 



To calculate the precession timescale, we have assumed a power-law profile for Bz, as in Lai 
(1998), such that: 



Bzir) = Bo,z[^J , (28) 

where x = correspond to the case of a constant Bz along the disk, while % = — 3 corresponds to 
the dipole case. 

Substituting our power-law parameterizations for S and Bz into equation (25) , we can calculate 
the relation between the observationally inferred precession period and Sq and Bq^z, expressed by 
the ratio taiupP^^BQ^/^o, as a function of the geometrical parameters a*, -Rout; * and %. 



3.4. Frame dragging and disk viscosity: the Bardeen-Petterson effect 

Frame dragging produced by a rotating compact body with angular momentum Jb causes 
precession of a particle if its orbital plane is inclined in relation to the equatorial plane of the 
rotating object. This is known as the Lense-Thirring effect (Lense & Thirring 1918). The precession 
angular velocity Olt produced by the frame dragging is given by (e.g., Wilkins 1972): 



2G Jb 



(29) 
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The presence of the Lense-Thirring effect in astrophysical systems with rotating neutron stars 

or Kerr black holes has been claimed in several works in the literature, as the physical mechanism 
behind the observed quasi-periodic oscillations (e.g., Cui, Zhang &; Chen 1998; Stella &; Vietri 1998; 
Markovic & Lamb 1998). 

The combined action of the Lense-Thirring effect and the internal viscosity of the accretion 
disk forces the alignment between the angular momenta of the Kerr black hole and the accretion 
disk. This is known as the Bardeen-Petterson effect (Bardeen k. Petterson 1975), and it tends to 
affect only the innermost part of the disk due to the short range of the Lense-Thirring effect, while 
the disk's outer part tends to remain in its original configuration. The transition radius between 
these two regimes is known as the Bardeen-Petterson radius, i?BP; and its location depends mainly 
on the physical properties of the accretion disk (Bardeen &; Petterson 1975; Kumar & Pringle 1985; 
Ivanov &: Illarianov 1997; Nelson & Papaloizou 2000). 

A rough estimate of i?BP can be obtained by comparing the time-scales for Lense-Thirring 
precession and warp transmission through the disk (e.g., Natarajan & Armitage 1999). If the 
transmission occurs diffusively, the Bardeen-Petterson radius can be obtained from: 



The diffusive regime in a Bardeen-Petterson disk has been explored by several authors (e.g., 
Kumar & Pringle 1985; Scheuer & Feiler 1996; Nelson &: Papaloizou 2000; Fragile & Anninos 2005; 
Lodato &; Pringle 2006), either using analytical calculations or numerical methods. 

A Bardeen-Petterson disk in a wave- like regime has also been studied in the literature (e.g., 
Ivanov & Illarianov 1997; Ogilvie 1999; Lubow, Ogilvie k, Pringle 2002); in such a situation, i?BP 
is given by: 



K 



BP 



cs(r = i^gp) 
^^LT(r = i?Sp) 



(31) 



The transition between the diffusive and wave-like regimes occurs approximately at a radius 
Rt ~ H/a (Papaloizou & Lin 1995). 

The time-scale for the black hole to align its angular momentum with that of the accretion 
disk was first estimated by Rees (1978). Scheuer & Feiler (1996) obtained an analytic solution to 
the equations that control the warp evolution in the case of a disk with constant surface density 
and used that to calculate the alignment time-scale. Natarajan k, Armitage (1999) generalized the 
results found by Scheuer k Feiler (1996) to a power-law viscosity. These studies suggest that the 
alignment time-scale can be estimated by: 
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Hiiign = Jb ( sin if, (32) 

where 99 is the angle between the angular momentum of the neutron star/black hole Jb and the 
direction perpendicular to the outer disk, and the time derivative of Jb is given by: 



/•Rout 

= -2iTsnnpj VLuY{r)LA{r)rdr. (33) 



Contrary to Scheuer & Fcilcr (1996), King ct al. (2005) showed that counter-alignment can 
occur when the spins are anti-parallel {(p > 7r/2) and the angular momentum on the disk is smaller 
than about twice that of the rotating compact object (Jd < 2Jb). The results obtained recently by 
Lodato Sz Pringle (2006) support the possible existence of counter-alignment. 

Other than in the special circumstances mentioned above, the Bardeen-Petterson effect forces 
the disk to align gradually with the rotating compact object. According to Scheuer & Feiler (1996), 
the time evolution of the alignment, and the precession period of the disk angular momentum around 
the spin axis of the rotating accreting object, -P^ec given by: 

ip{t) = ipoe-^^/^-^'s- (34) 

P^ii) = Poe-^*/^-- (35) 

where At = t — to, and ipo and Pq are, respectively, the inclination angle and precession period at 
time to , when the action of the Bardeen-Petterson torques over the accretion disk started (to < is 
measured in the past from the present time). Scheuer & Feiler (1996) found that the timescales for 
precession and realignment are identical, implying that Pq = T'aiign, which will also be used here. 

4. Source sample 

To analyze separately the consequences of each of the precession mechanisms discussed in the 
last section, we selected eight X-ray binaries and four AGNs that present signatures of warping 
and/or precession in their accretion disks. In this section we introduce these systems and provide 
their basic parameters (the ones that will be used in our calculations). 

X-ray binaries, characterized as semi-detached systems with one component filling its critical 
Roche lobe, can be roughly divided into two categories: high- and low-mass X-ray binaries (HMXBs 
and LMXBs respectively). 
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HMXBs are composed of an accreting compact source (a black hole or a neutron star) and 

either an OB supergiant or a Be star (van Paradijs 1983), which implies that such systems are 
short-lived, with ages of less than about 2 x 10^ yr (van den Hcuvcl 1994). They are spatially 
distributed in the galactic plane, being associated with young stellar populations. 

By contrast, LMXBs contain a neutron star accreting from a low-mass star (< 2Mq). They are 
typically associated with an older stellar population, with ages of about one Gyr (van den Heuvel 
1994). 

Here, we will use eight X-ray binaries with regular (or quasi-regular) precession periods inferred 
basically from the optical and X-ray continuum variability, as well as from variability in the line 
intensity and velocity in some cases. Our sample is composed of two LMXBs (HerX-1 and CygX- 
2) and six HMXBs, four of which have an accreting neutron star (LMCX-4, CenX-3, SMCX-1 
and 4U 1907-1-09) while the remaining sources are black hole systems (SS433 and CygX-1). It is 
important to note that SS433 is one of the best studied microquasars in the literature, with a 
prominent precessing jet/counterjet (e.g., Hjellming & Johnston 1981; Margon 1984; Blundell et 
al. 2001; Blundell k. Bowler 2004). The relevant parameters of the sources, for the present work, 
are listed in Table 1. In particular, the values of the magnetic field strength listed in Table 1 refer 
basically to the surface magnetic fields of neutron stars, while for the two black hole systems we 
used values derived from dynamo models for the disk magnetic field (Pudritz &; Fahlman 1982; 
Rose 1995). 

Some AGNs also exhibit signatures of warping/precession in their accretion disks. The signa- 
tures are in the form variability of double-peaked Balmer lines and of the continuum emission, or 
distortions in the jet morphology (variations in the jet orientation and jet velocity) (e.g., Abraham 
2000; Storchi-Bergmann et al. 2003; Caproni & Abraham 2004b). We have selected four AGNs with 
these characteristics: the Seyfert 1 galaxy 3C 120, the Seyfert 2 galaxy NGC 1068, the broad-line 
radio galaxy Arp 102B, and the BL Lac OJ 287. Their parameters are listed in Table 2. In relation 
to their magnetic field strengths, since there is evidence in some cases that the strength of the 
magnetic field is not in equipartition (e.g. Coker & Melia 2000), we have chosen the field strength 
to be equal to that which has been observationally determined in a few AGN (~ 10^ G; e.g.. Field 
& Rogers 1993; Lobanov 1998). 

5. Analysis 

We have analyzed which of the four mechanisms can be driving the observed precession in 
the twelve sources of our sample. Below we describe the results of this analysis, discussing each 
mechanism separately. 
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5.1. Tidal torques in a binary system 

In order to examine the possibility that tidal torques induced by the companion affect the 
accretion disk, we have plotted in the left panel of Figure 1 the quantity cos 0Pprec/-Porb predicted 
by equation (14) as a function of q (thick line); the values for all the X-ray binaries in our sample 
are also displayed (star symbols). The values of 9 for HerX-1, LMCX-4, SS433 and CygX-1 were 
taken from the literature (Heemskerk & van Paradijs 1989; Scott et al. 2000; Stirling et al. 2002; 
Romero et al. 2002), while 6 = 10° was assumed for the other X-ray binaries. 

The gray regions in Figure 1 show the part of parameter space where tidal torques from a 
companion can drive precession. Clearly, the observed ratios between the precession and orbital 
periods for CenX-3 and CygX-2 reside below of this region. Therefore, the origin of their disk 
precession cannot be attributed to the tidal torques of the secondary. In the case of CygX-1, 
Lachowicz et al. (2006) have recently claimed that its observed precession is prograde, which would 
also rule out the tidal torque scenario. For the other five X-ray binaries, precession may be indiiccd 
by the torques of the companion star. Although Katz (1973) and Larwood (1998) have already 
proposed that tidal torques could lead to the observed precession in a few of these sources, tidally- 
induced precession has not been previously suggested (to the best of our knowledge) for 4U 1907+09. 

In the case of the AGN sources, there is no straightforward observation that supports the 
existence of a supermassive binary black hole system in any of their nuclear regions. Nevertheless, 
the continuum variability and the anomalous jet kinematics in O J 287 and 3C 120 have been used 
to put some constraints on the physical parameters of a possible binary system (Sillanpaa et al. 
1988; Abraham 2000; Caproni & Abraham 2004b). In the right-hand panel of Figure 1, we show 
the upper limit for function of q for our extragalactic sources, considering a 

timescale for losses due to gravitational radiation tgw (Shapiro & Tcukolsky 1983) of 1000 yr. 
Although this value has been chosen arbitrarily, it guarantees that no significant changes in the 
orbit of the secondary occur on such a timescale, so that the observed precession periods do not 
vary substantially. Note that increasing (decreasing) tqw by a factor of 10 results in a decrease 
(increase) in the upper limit for cos 0Pprec/-Porb only by a factor of ~2.4. For OJ287, 3C 120 and 
NGC 1068, 9 was taken from the literature (Abraham 2000; Caproni Sz Abraham 2004b; Caproni, 
Abraham Sz Mosquera Cuesta 2006), while for Arp 102B, we assumed 9 = 10°. 

By contrast with the galactic binaries, we do not expect g > 1 in supermassive binary black 
hole systems, since the expectation is generally that the accretion disk responsible for the AGN 
activity is associated with the more massive black hole. Except for Arp 102B, it is always possible to 
find a value of q <1 for which cos ^Pprec/-Porb is within the gray area. In addition, the separations 
between the putative black holes in NGC 1068 obtained from our calculations ^ are always smaller 
than the dimensions of the mascr disk. Consequently, we would expect a more complex disk 
morphology than that suggested by the observations. Thus, tidal forces of an orbiting black hole 



^Typically 100-10'* ilg, depending on the timescales for losses due to gravitational radiation. 
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can (in principle at least) induce the inferred precession rate in 3C 120 and OJ 287 (and perhaps 
NGC 1068). 

5.2. Irradiation-driven torques 

In Figures 2 and 3 wc present the precession period induced by the radiation torques as 
a function of the critical radius at which the disk becomes irradiation-warping unstable, for each 
source of our sample, using equations (15) and (16). Following Maloney, Begelman & Nowak (1998), 
the calculations were performed for power-law disks, using s = —3/2, —1, and 3/2. The lower and 
upper limits for each solution refer to the respective lower and upper limits of Sq calculated from 
equations (9), (11) and (19), as discussed previously. 

The critical radius, as well as the precession period, depend on the ratio between rj and 
e, which is usually unknown. We have assumed 0.0377 < e < 0.42 ^, while 1 < rj < f{a), 
where /(a) = 2(1 + 7a^)/[a^(4 + a^)] (Ogilvie 1999). In the case of ry = 1 (open symbols in 
Figures 2 and 3), 2.38 < r//e < 26.49, implying 223.80 < Rcr/Rg < 2.77 x 10*^ for s = -3/2 and 
1.34 X 10^ < Rcr/Rg < 1.66 X 10^ for s = 3/2. For r? / 1, ??/e, and consequently P^f^^ and i?cr, 
also depend on the value of a; we have chosen a = 0.1 for the X-ray binaries and a = 0.01 for the 
AGNs, except for NGC 1068 for which 0.001 < a < 0.012 (Caproni, Abraham & Mosquera Cuesta 
2006). These values of a imply that the critical radii, for a fixed e, are systematically larger in the 
case of 77 7^ 1. 

For the X-ray binaries, in the case of an isothermal disk (s = —3/2) with 77 = 1 (open circles), 
the predicted precession period intervals cross the horizontal line, which represents the observed 
precession period. This indicates that the radiation instability can be responsible for precession in 
those X-ray binaries (for this particular accretion disk model). However, the situation changes when 
we take into account other combinations among s, rj and e; even though the predictions assuming 
s = —3/2 and r) 1 are compatible with the observed precession in most cases, power-law disks 
with s = 3/2 generally fail to reproduce the observations, providing incompatible precession periods 
and/or critical radii larger than the outer disk radii. 

As pointed out by Maloney, Begelman &; Nowak (1998), disk precession induced by radiation 
must be prograde in the absence of external torques, which would be in contradiction with the 
observations of SS433, HerX-1 and LMCX-4. However, as shown by Maloney & Begelman (1997) 
and corroborated by Maloney, Begelman & Nowak (1998), the quadrupole torque from a companion 
star might allow for the existence of prograde and retrograde precession modes. Thus, retrograde 
precession due to the irradiation mechanism is possible only if there is some additional torque acting 
upon the accretion disk. 



^The upper and lower limits correspond to the accretion efficiencies from the maximum prograde and retrograde 
spinning black holes, respectively. 
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In the case of AGNs, the situation is quite different. For r) ^ 1, the critical radius is larger 
than the outer disk radius for all AGNs in our sample, suggesting that in these cases the disk is 
stable against radiation torques. 

Thus, the irradiation instability can be responsible for the precession in our sample of X-ray 
binaries, in agreement with Ogilvie k, Dubus (2001), but it is not favored in the case of the AGNs 
we have considered. 

5.3. Magnetically-driven torques 

To analyze the precession induced by magnetic torques, we created a two-dimensional grid of 
the quantity tan tpP^ei ^/^o^ fixing s and x, as a function of a* and Rout for each source in our 
sample. 

These grids arc plotted in Figures (4) and (5) for the X-ray binaries and in Figure (6) for the 
AGNs assuming power-law disks with s = —2, —1 and and power- law magnetic field configurations 
with X = 0, —1, —2 and -3. The gray area represents the allowed range for P^^B^ 2/^0 (assuming 
tamp = 1) using the available constraints obtained from the observations. As in the case of tidal 
torques in binary systems, acceptable solutions must also obey -Rprcc < -Rout- The dimensionless 
spin parameter influences only the location of the inner radius of the accretion disk, while the 
estimated extreme values for Sq are responsible for the limits of the gray region in the panels of 
the Figures (4)- (6). 

We can see in Figures (4)- (6) that not all magnetic field configurations provide solutions 
compatible with the observations. Nevertheless, it is possible to find combinations of s and x 
that reproduce the observed precession periods, which means that magnetically-driven instabilities 
cannot be excluded as a potential mechanism for the observed precession. 

At least in the case of the X-ray binaries with an accreting neutron star, such as HerX-1, 
LMCX-4, CenX-3, CygX-2, SMCX-1 and 4U 1907-1-09, the magnetic field configuration is dipolar 
to a good approximation (e.g., Aly 1980). We have not found any possibility that disk precession in 
those systems can be associated with magnetic field configurations with x ^ "Ij except for CygX-2 
for which acceptable solutions are obtained with x = — 1 and x = only if s = —2. For the black 
hole X-ray binaries, CygX-1 has similar results to those of CygX-2, while for the microquasar 
SS 433 lower values of x are strongly favored. 

In contrast to the neutron star X-ray binaries, there is no preferred magnetic field configuration 
in the case of AGNs. Our results show that the range of the allowed solutions is systematically 

narrower than that obtained for the X-ray binaries. In addition, it seems to have two different 
regimes in our AGN sample: if the accretion disk surface density decreases faster with radial 
distance, with s < — 1, the z-component of the magnetic field must be radially constant in order to 
reproduce the observed precession periods; otherwise, it is necessary for the surface density not to 
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vary radially (s = 0) if the magnetic field decreases along the disk (x < — 1). 



5.4. Bardeen-Petterson efTect 

In order to study the consequences of the Bardeen-Petterson effect in our sample of sources, 
we have separated the black hole systems from those with accreting neutron stars. In the former 
case, we have calculated the Bardeen-Petterson radius and the alignment timcscale for six values 
of the black hole spin (a^, = ±0.1, ±0.5 and ±1) and for s = —2, —1 and 0, as shown in Figure 7. 
As in the previous sections, calculations were performed for rj = f{a) (a = 0.1 and 0.01 for X-ray 
binaries and AGNs, respectively) and = 1, as well as for the estimated lower and upper limits on 
the accretion disk surface density. 

Once s and r] are fixed, there are generally two values of Rsp that are determined by the lower 
and upper limits of Sq. In some cases the lower limit for Sq leads to Rsp < Rms, and because of 
that only the upper limit is shown in the figure. All sources have their Bardeen-Petterson radius 
inside of the outer edge of the disk, independently of the specific values for s, rj and Sq, which 
indicates that the Bardeen-Petterson mechanism is at least applicable. 

We also plot in Figure 7 the time interval necessary for each system to reach the observationally 
inferred precession period (using equation 35). As in the case of the alignment timescales, the plots 
were truncated at 13.7 Gyr, the age of the Universe as inferred from the WMAP results (Bennett 
et al. 2003). Unfortunately, only in the case of SS433 we can constrain better the allowed model 
parameters using its estimated age (between 2 x 10^ and 2 x 10^ yr; Zealey, Dopita, &; Malin 
(1980)). This age constraint clearly favors model disks with rj ^ 1 at least for |a^<| < 0.1. The same 
conclusions are valid for CygX-1 if we take into account the upper limit for the age of HMXBs 
(van den Heuvel 1994). 

In the case of NGC 1068, we also have included the solutions obtained from a power-law 
accretion disk with s = —1.05 (Hure 2002), corresponding to model A shown in Figure 2 of Caproni, 
Abraham & Mosquera Cuesta (2006). Indeed, these authors showed that the Bardeen-Petterson 
effect can reproduce the disk configuration suggested by the maser observations (Gallimore et al. 
2004), as well as the general shape of the parsec and kiloparsec radio jet. 

For the neutron star systems in our sample, the spin period of the accreting neutron stars, Pg, 
are known, and are listed in Table 3. From the definition of the angular momentum parameter, we 
can obtain (e.g., Stella & Vietri 1998): 

27rc Ip Us 
= iWp Mp 

where Ip and are respectively the moment of inertia and the spin frequency of the neutron star. 
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If we write Ip = 10^^/45 g cm^ and Mp = Mq M© typically, 0.5 < /45/Mo < 2, considering 
rotating neutron star models for different equations of state and masses. In Table 3 we give the 
allowed ranges of a* for our neutron star systems using their respective orbital periods and the 
lower and upper limits for /45/Mo. 

In Figure 8, we plot the results for the neutron star system displaying, as in the Figure 7, the 
Bardeen-Petterson radius, the alignment timescale and the time interval necessary for each system 
reach the observationally inferred precession period for the maximum absolute value of a* listed in 
Table 3. As in the case of the Kerr black holes, all sources have their respective Bardeen-Petterson 
radii inside the outer radius of the accretion disk. All spinning neutron stars in our sample reach the 
observed precession periods in a timescale shorter than 20 Myr. In the case of Her X-1, whose age 
has been estimated previously (Verbunt, Wijers & Burm 1990), any allowed combination of model 
parameters leads to Her X-1 reaching its observed precession period via the Bardeen-Petterson 
effect. The situation is similar, even though somewhat more restrictive, for CygX-2. 



6. Conclusions 

In this work, we have selected eight X-ray binaries and four AGNs that present signatures 
of precession in their accretion disks in order to examine the compatibility between their preces- 
sion periods and the predictions from four distinct warping/precession physical mechanisms: tidal 
torques from a companion in a binary system, irradiation- and magnetically-driven instabilities, 
and the Bardeen-Petterson relativistic effect. 

We have assumed a power-law surface density distribution for the accretion disks, constraining 
their physical parameters from observational data available in the literature. 

For the X-ray binaries in our sample, we found that tidal torques from a companion in a binary 
system provide precession timescales compatible with those inferred in SS433, Her X-1, LMCX-4 
and SMC X-1, as indeed has been previously suggested in the literature (Katz 1973; Larwood 1998). 
In addition, we showed that tidal torques can also drive precession in 4U 1907-1-09, which (as far as 
we know) had not been proposed before this work. The mechanism can be ruled out for CenX-3 
and CygX-2. In the case of our AGN sample (assuming that they contain binary black holes), we 
have shown that tidal torques cannot produce the disk precession inferred in Arp 102B. 

Although the irradiation-driven instability usually provides precession timescales that agree 
with those observed in our X-ray binary sample, it is unsuccessful in reproducing the reported pre- 
cession periods for the four AGNs considered here. The critical radii at which those AGN accretion 
disks become unstable against radiation torques were found to be larger than the expected disks' 
outer radii (for / 1). Since the irradiation-driven instability produces prograde disk precession, 
the retrograde precession observed in SS433, Her X-1 and LMCX-4 may only be generated in the 
presence of some external torques. 
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We have shown that torques due to a misahgned magnetic field (with respect to the perpen- 
dicular direction to an accretion disk) can also reproduce the observed precession periods of our 
sample. For the accreting neutron stars, a dipolar configuration must be a good approximation 
to their magnetic fields; in such a case (or for x < — 1), the magnetic-driven torquing can induce 
precession in the inner parts of the accretion disks at the rate inferred observationally. For the black 
hole X-ray binaries, magnetic configurations with x ^ ~1 favored. Although the geometry of 
the magnetic field is not usually constrained by observations in AGNs, our results indicate the 
existence of two distinct regimes: for accretion disks with s < —1, the z-component of the magnetic 
field must be radially constant in order to reproduce the observed precession periods; otherwise, it is 
necessary to have a constant surface density (s = 0) with a magnetic field weakening radially along 
the disk (x < — !)• In addition, if the accretion disk of CygX-1 is actually precessing progradely 
(Lachowicz et al. 2006), we can rule out magnetically-driven torques as the cause of its precession. 

The Bardeen-Petterson effect produces precession timescales compatible with those observed 
in all sources of our sample, considering that the alignment of the accretion disk evolves on a 
similar timescale (as in Scheuer & Feiler 1996). The timescale for reaching the observed precession 

rate is usually shorter than the estimated lifetime of the sources. Considering the accretion disk 
parameters used in this work, we have found that all sources in our sample have Bardeen-Petterson 
radii smaller than the disk's outer radius. 

NGC 1068, for which maser observations have been used to infer the physical properties of the 
accretion disk, may provide the strongest candidate for the Bardeen-Petterson effect. In this case, 
we have shown that the irradiation-driven instability is incompatible with the observations. Tidal 
torques are also highly unlikely (even if a binary black hole were to be present in this system) 
since it would most likely result in other changes to the disk morphology, which are not observed. 
Magnetically-driven torques would work only if a rather contrived field configuration exists {Bz 
constant along the precessing part of the disk) . Our conclusions give further support to the analysis 
of Caproni, Abraham & Mosqucra Cuesta (2006), who showed that the Bardeen-Petterson effect is 
consistent with the warping morphology suggested by Gallimore et al. (2004). 

The Bardeen-Petterson precession is an important general-relativistic effect, and its unam- 
biguous detection would be of great value. This work has shown, however, that given the many 
observational uncertainties that still exist (both in X-ray binaries and in AGNs), it is extremely 
difficult to rule out definitively other precession mechanisms. 

It is important to emphasize that although we have analyzed separately the four precession 
mechanisms in this work, they are not mutually exclusive. This means that more than one mecha- 
nism might be operating in a given system. 

A.C. acknowledges financial support from FAPESP (Procs. 03/13882-0 and 05/56266-2), as 
well as the hospitality of the Space Telescope Science Institute, where this work was carried out. 
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Table 1. Physical parameters of the X-ray binary sources. 



Source 


Mp (Mq) 


Ref. 




Porh (d) 


Ref. 


-fprcc (d) 


Ref. 


Type"! 


Ref. 


B (G) 


Ref. 


Lbol (erg s 


Ref. 


SS433 


11 ± 3 


1 


1.73 


13.08 


8 


162.5 


15 


R 


23,24 


10^ 


28 


39.74 


36 


HerX-1 


1.50 ± 0.30 


2 


1.53 


1.7 


9 


34.88 


16 


R 


25 


3.5x10^2 


29 


37.58 


2 


LMCX-4 


1.38 ± 0.25 


3 


10.6 


1.408 


3 


30.275 


17 


R 


26 


10" 


30 


38.19 


37 


CcnX-3 


1.21 ± 0.21 


4 


16.9 


2.087 


10 


140=^ 


18 


7 




3.5x10^2 


31 


37.88 


38 


CygX-1 


10.1 


5 


1.76 


5.67 


11 


142 


19 


P 


27 




32 


37.80 


39 


CygX-2 


1.78 ± 0.23 


6 


0.34 


9.844 


12 


77.7 


20 


? 




4.5x10*= 


33 


37.65 


40 


SMCX-1 


1.60 ± 0.10 


7 


10.75 


3.893 


13 


60 


21 


? 




7.6x10^2 


34 


38.30 


41 


4U 1907+09 


1.4*^ 




21.4 


8.375 


14 


41.6 


22 


? 




2.1x10^2 


35 


37.728 


42 



References. — (1) Gies, Huang & McSwain 2002; (2) Reynolds et al. 1997; (3) Levine et al. 1991; (4) Ash et al. 1999; (5) Herrero et al. 1995; (6) Orosz 
& Kuulkers 1999; (7) Reynolds et al. 1993; (8) Stewart et al. 1987; (9) Tananbaum et al. 1972; (10) Paul, Raichur & Mukherjee 2005; (11) Priedhorsky, 
Brandt & Lund 1995; (12) Casarcs, Charles & Kuulkers 1998; (13) Schrcicr ct al. 1972; (14) in't Zand ct al. 1998; (15) Stirling ct al. 2002; (16) Giacconi 
ct al. 1973; (17) Tsygankov & Lutovinov 2005; (18) Priedhorsky & Terrell 1983; (19) Brocksopp ct al. 1999; (20) Wijnaiids, Kuulkers & Smalc 1996; (21) 
Wojdowski et al. 1998; (22) Priedhorsky & Terrell 1984; (23) Lcibowitz 1984: (24) Brinkmann, Kawai & Matsuoka 1989; (25) Gcrcnd & Boynton 1976; 
(26) Heemskerk & van Paradijs 1989; (27) Lachowicz et al. 2006; (28) Rose 1995; (29) dal Fiume et al. 1998; (30) La Barbera et al. 2001; (31) Burderi 
et al. 2000; (32) Pudritz & Fahlman 1982; (33) Campana 2000; (34) Kahabka & Li 1999; (35) Coburn et al. 2002; (36) Okuda 2002; (37) Vrtilek et al. 
1997; (38) Lanzafame & Belvedere 1998; (39) Zdziarski et al. 2002; (40) Vrtilek et al. 1990; (41) Wojdowski, Clark & Kallman 2000; (42) Mukerjee et al. 
2001. 

"Mass value assumed in this work; 

''Quoted errors obtained from error propagation; 

"^Mean value; 

''Type of precession (in relation to the angular momentum of the accretion disk: prograde (P), retrograde (R). The symbol "?" indicates that there is 
no available information; 

''Upper limit; 

'Base-10 logarithm of the bolometric luminosity; 
^Flare state (upper limit). 
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Table 2. Physical parameters of the AGNs. 



Source 


Mtot (10* Mo) 


Ref. 


Ppvec (yr)'' 


Ref. 


Type-^ 


B (G)d 


Lhol (erg s 


Ref. 


NGC 1068 


0.120 


1 


(2.08-35.2) xlO^'' 


5 


7 


10'' 


44.84' 


9 


3C 120 


0.556 


2 


11.9 


6 


? 


10" 


45.34 


10 


OJ287 


6.170 


3 


8.88 


7 


7 


10" 


44.448 


3 


Arp 102B 


1.380 


4 


2.15 


8 


7 


10" 


43.41 


4 



References. — (1) Hure 2002; (2) Peterson et al. 2004; (3) Wang, Luo & Ho 2004; (4) Wu & Liu 2004; (5) 
Caproni, Abraham & Mosquera Cuesta 2006; (6) Caproni & Abraham 2004b; (7) Abraham 2000; (8) Newman et 
al. 1997; (9) Gallimore et al. 2004; (10) Woo & Urry 2002. 

^Values measured at the present time and in the source's reference frame; 

''Obtained from the jet kinematics and the Bardeen-Petterson model (Caproni, Abraham & Mosquera Cuesta 
2006); 

■^Same nomenclature of the Table 1; 

''Adopted value for the magnetic field at Rg; 
''Basc-10 logarithm of the bolometric luminosity; 
'This value is a lower limit; 

^Assuming also that the bolometric luminosity is about ten times larger than the luminosity of the broad line 
region (Netzer 1990). 



Table 3. Rotation period and spin parameter for the neutron stars in our X-ray binaries sample. 



Source 


(s) 


Ref. 






HerX-1 


1.2377697 ± 0.0000003 


1 


(1.92-7.69) X 10" 


-4 


LMCX-4 


13.509 ± 0.002 


2 


(1.91-7.66)xl0- 


-5 


CenX-3 


4.834477 ± 0.000007 


3 


(0.61-2.44) X 10" 


-4 


CygX-2 


0.00606 ± 0.00024 


4 


(0.33-1.32) X 10^ 


^1 


SMCX-1 


0.706707 ± 0.000001 


5 


(0.32-1.26) X 10" 


-3 


4U 1907+09 


440.5738 ± 0.0002 


6 


(0.58-2.31) X 10" 


-6 



References. — (1) Oosterbroek et al. 2001; (2) Vrtilek et al. 1997; 
(3) van der Klis, Bonnet-Bidaud & Robba 1980; (4) Focke 1996; (5) 
Kahabka & Li 1999; (6) Baykal et al. 2001. 

*'The lower and upper limits refer respectively to li^/Mo = 0.5 and 
2 (Stella & Vietri 1998). 



-28- 




Fig. 1. — Precession due to tidal torques induced by the companion in binary systems. Left panel: 
Ratio between the precession and orbital periods (multiphed by the cosine of the incUnation angle 
of the orbital plane in relation to the disk plane) as a function of the mass ratio between the 
secondary and primary objects. The solid line is the theoretical prediction considering that the 
entire disk with an outer radius of 0.88i?L precesses rigidly. Data from our sample of X-ray binaries 
are displayed by the stars. Sources located in the gray area can have their disk precession driven 
by the companion's torque. Right panel: Upper limit for the AGNs of our sample considering a 
timescale for losses due to gravitational radiation equal to 1000 yr. The dashed vertical line marks 
q = l. 
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Fig. 2. — Precession due to the radiation torques in our sample of X-ray binaries. Circles, triangles, 
squares and diamonds correspond respectively to s = —3/2, —1, and 3/2. Filled symbols represent 
the case 77 = /(a) for a = 0.1, while the open ones 77 = 1. The solid horizontal line in each 
panel refers to the observationally inferred precession period, while the vertical solid line mark the 
location of the outer disk radius. The allowed range for the precession period is represented by the 
connecting lines, whose extremes are given by the lower and upper limits of the surface density of 
the accretion disk at the gravitational radius. Model predictions beyond the outer disk radius (not 
physically acceptable) are plotted as dashed lines. 
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Fig. 3. — Precession due to the radiation torques applied to our sample of AGNs. Same nomen- 
clature as in Figure 2 is assumed here. In the case of rj = f{a), we have chosen 0.001 < a < 0.012 
for NGC 1068 (Caproni, Abraham & Mosquera Cuesta 2006) and a = 0.01 for the other sources. 
The crossed gray squares in the panel of NGC 1068 refer to the modeling of the accretion disk by 
Hure (2002) based on maser data. The solid horizontal line refers to the observationally inferred 
precession period. The vertical lines refer to the disks' outer radii, 1.1 pc for NGC 1068 (inferred 
from interferometric maser observations Gallimore et al. 2004) and 10^i?g for the other sources 
(upper limit based on Collin- Souffr in & Dumont 1990). 
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Fig. 4. — Magnetically-driven instability. The logarithm of tan(^Pj^|f Z?q ^/Sq is calculated as a 
function of the dimensionless spin parameter a* and the outer radius of the warped/precessing part 
of the disk i?prcc (normalized by the radius of the marginally stable orbit) for SS433, HcrX-1, 
LMCX-4 and CenX-3. P^cc is given in years, while i?o,z and Sq are in CGS units. The gray 
region in the panels correspond to the allowed range constrained by the observations considering 
tan(p = 1. Calculations were performed (see text) for s = —2, -1 and (from the left to the right 
side) and for x = 0, -1, -2 and -3 (from the top to the bottom panel). 
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Fig. 5. — Magnetically-driven instability for the galactic sources CygX-1, CygX-2, SMCX-1 and 
4U 1907-1-09. Same nomenclature as in Figure 4 is adopted here. 
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Fig. 6. — Magnetically-driven instability for the extragalactic sources NGC 1068, 3C 120, OJ287 
and Arpl02B. Same nomenclature as in Figure 4 is adopted here. The thick black line represents 
the solution for NGC 1068 found from the accretion disk model parameters given by Hure (2002) 
(see also Caproni, Abraham &; Mosquera Cuesta 2006). 
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Fig. 7. — Bardccn-Pcttcrson effect in accreting black holes. We show the Bardeen-Petterson radius, 
the timescale for the alignment between the angular momenta of the black hole and of the accretion 
disk, and the time interval for the system to reach the observationally inferred precession period. 
Circles, triangles and squares represent respectively the solutions for s = —2, -1 and 0. Filled 
symbols correspond to = /(a) {a = 0.1, for SS433 and CygX-1, and a = 0.01 for the AGN 
sources) while the open ones are for rj = I. Stars in the panels of NGC 1068 show the results from 
the disk model parameters given by Hure (2002) (see also Caproni, Abraham & Mosquera Cuesta 
2006). The horizontal lines in the panels on Rbp show the values of the disk outer radius. The 
two parallel lines in the Af-plot for SS 433 are the estimated age of SS 433 (between 0.02 and 0.2 
Myr; Zealey, Dopita, &; Malin 1980), while the dashed line seen in the At-plot for CygX-1 is the 
approximated upper limit for the ages of HMXBs (van den Heuvel 1994). The upper limit for age of 
NGC 1068 is 3.5 x 10^ yr (Wilson k Ulvestad 1987). All plots referring to timescales are truncated 
at 13.7 Gyr. 
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Fig. 8. — Bardeen-Petterson effect in accreting neutron star in X-ray binaries. Same nomenclature 
as in Figure 7 is adopted liere. The lower limit for tfie age of HerX-1 is ~ 6 x 10^ yr (Verbunt, 
Wijers & Burm 1990), while the estimated age of CygX-2 is roughly 3 Myr (Kolb et al. 2000). 



